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Order a improved renormalization constants 
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We impose the axial and vector Ward identities on local fermion bilinear operators in the Sheikholesalami- 
Wohlert discretization of fermions. From this we obtain all the coefhcients needed to improve the theory at 0{a), 
as well as the scale and scheme independent renormalization constants, Za, Zy and Zs/Zp. 
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1. Introduction 

In Reference [0, we studied the improvement 
of the Wilson-Dirac theory by removing all lattice 
artifacts linear in the lattice spacing a in on-shcU 
matrix elements. In general, such an improve- 
ment for the dimension 3 fermion bilinear oper- 
ators requires one to tune the coefficient of the 
clover operator in the action, determine the coef- 
ficient of an extra term each in the vector, axial- 
vector, and tensor currents, and obtain the mass 
dependence of all the renormalization constants. 
We showed that all these improvement constants, 
as well as the renormalization constants for the 
axial and vector currents, can be determined by 
imposing the axial and vector Ward identities. 
The remaining renormalization constants, those 
for the scalar, pseudoscalar and tensor operators, 
are scheme dependent and cannot be determined 
by this method. The scheme dependence is, how- 
ever, the same for both the scalar and the pseu- 
doscalar operators, and the corresponding ratio of 
renormalization constants can again be obtained. 

To test the efficacy of this method, we had eval- 
uated the renormalization constants at /3 = 6.0 
for the perturbative (tadpole improved) value of 
the clover coefficient, and found that these dif- 
fered from the values determined by the ALPHA 
collaboration ||] using a non-perturbatively tuned 
coefficient. It was not clear whether the differ- 
ences arose from different O(a^) errors in the two 
calculation, or whether the difference in the clover 
coefficient was responsible for the variation ob- 



served. Here we report results of carrying out our 
procedure after changing the clover coefficient to 
that used by the ALPHA collaboration, and also 
at /3 = 6.2, where the ©(a^) errors are expected 
to be much smaller. 

For the details of the notation used, the numer- 
ical technique, a study of the consistency when 
the same constant is determined in multiple ways, 
and the choice of the method which determines 
each quantity best, we refer the reader to Ref- 
erence To avoid confusion, we repeat here 
that the coefficients hx we defined differ from the 
coefficients hx used by earlier authors. In partic- 
ular, at the level of 0{a) improvement, one has 
bx = {Z\Zl/Z%)hx. 

2. Lattice parameters and Results 

The results of our calculation done at /3 = 6.0 
and 6.2 are presented in Table |l| and compared to 
the previous estimates by the ALPHA collabora- 
tion dl and perturbation theory. In those cases 
where the coefficients hx and hx differ by more 
than the estimated error in our determination, we 
quote both for easier comparison with previous 
results. 

The (3 = 6.0 calculations were done on 16'^ x 
48 lattices, with a sample of 83 configurations at 
csw — 1.4755, and 125 configurations at csw = 
1.769. At /3 = 6.2, we used 61 configurations 
with a lattice size of 24^ x 64. The region of axial 
rotation used for the Ward identities comprise the 
time-slices 4 — 18 at /? = 6.0 and 6 — 25 at /3 = 
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Table 1 



Renormalization constants and improvement coefficients at (3 — 6.0 and /3 = 6.2. 
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+0.745(1)(1) 
+0 76f2)fl) 
+0.77(4)(1) 


+0.764(1)(1) 

+0.811(9)(12) 

+0.86(2)(G) 


G.78G9(6) 
G. 7906(94) 
N.A. 


+G.81G(5) 
+G. 829(3) 
+G.949(3) 


+0.784(00)(00) 
+0.816(08)(05) 
+0.878(15)(06) 


+0.7922(4)(9) 

+0.807(8)(2) 

N.A. 


+0.821(5) 
+0.839(2) 
+0.953(2) 


CA 
CV 
CT 


-G.005(15)(20) 
-0.66(27)(02) 


-O.G2(2)(2) 
-O.G2(15)(G7) 


-O.G83(5) 
-0.32(7) 
N A 


-0.013(0) 
-0.028(1) 
-1-0 r\9r\(r\\ 


-0.033(04)(03) 
-0.10(14)(06) 
-1-0 of;9(''?^vns'i 


-0.038(04) 
-0.21(7) 
M A 


-0.012(1) 
-0.026(7) 
01 oi(A\ 


bv 

b A - by 

bp - bs 
bp - Ia 


+1.57(2)(1) 
+1.62(3)(1) 
-G.56(9)(8) 
-0.09(9)(3) 
-0.20(9)(7) 


+ 1.62(2)(2) 
+ 1.68(3)(3) 
-0.26(6)(9) 
-O.G2(6)(G) 
-O.G9(8)(6) 


N.A. 
+ 1.54(2) 
N.A. 
N.A. 
N.A. 


+1.106(11) 

+1.273(16) 

-0.002(0) 

-0.066(4) 

+0.002(4) 


+1.48(1)(1) 
+1.54(1)(0) 
-0.11(4)(3) 
-0.06(3)(1) 
-0.07(2)(3) 


N.A. 
+ 1.41(2) 
N.A. 
N.A. 
N.A. 


+ 1.099(10) 
+ 1.254(14) 
-0.002(0) 
-0.062(4) 
+0.002(0) 


Ia - bp 
+ bs/2 


+0.48(1)(2) 


+0.69(25)(1) 


N.A. 


+0.584(7) 


+0.50(5)(0) 


N.A. 


+0.579(6) 


bA 
bp 
bs 


+1.01(09)(09) 
+0.81(14)(16) 
+0.90(17)(13) 


+1.36(G6)(1G) 
+ 1.28(G9)(15) 
+1.29(1G)(12) 


N.A. 
N.A. 
N.A. 


+1.104(11) 
+1.105(11) 
+1.172(30) 


+1.37(4)(2) 
+1.30(4)(5) 
+1.36(5)(4) 


N.A. 
N.A. 
N.A. 


+ 1.097(9) 
+ 1.099(10) 
+ 1.161(26) 



6.2. Details of the quark masses and the lattice 
discretizations employed in the calculation will be 
presented elsewhere 

3. Dependence on (3 and csw 

The determination of ca is illustrative of the 
improvement we notice as we go to weaker cou- 
plings. This quantity is determined by requiring 
that the ratio 

T.s{dM^ + acAd^Pf^\S,t)j'^^^HQ)) 

be independent of the time t at which it is evalu- 
ated, up to errors of O(a^). Because this flatness 
criterion is automatically satisfied if the correla- 
tors are saturated by a single state, the determi- 
nation is very sensitive to the small t region. On 
the other hand, at very small t, the O(a^) errors 
dominate. We find that the region well fit by a 
constant starts at smaller values oi t e.t (3 — 6.2, 
and that this results in smaller statistical errors 
on the determination of ca- As ca feeds into the 



3-pt Ward identity calculation, the statistical er- 
rors are typically smaller at the weaker coupling. 
In addition, the sensitivity to how the contin- 
uum derivatives are discretized on the lattice is 
reduced, leading to smaller systematic errors. 

In contrast, the proper choice of csw at /3 = 6.0 
has a much smaller effect on the quality of the sig- 
nal. We find that changing csw from 1.4755 to 
1.769 has little effect on the statistical and sys- 
tematic errors but brings Z^^ and cx closer to 
their perturbative values. In addition, the differ- 
ences between the various bx, which are almost 
zero in perturbation theory, also decrease. 

4. Comparison with previous results 

Different calculations of the non-perturbative 
renormalization constants are expected to differ 
because of residual O(a^) errors in the theory. 
The magnitude of these effects are expected to 
be 0{AQfjjja'^) in Z^ and 0{Aqcdci) in cx and 
bx- Numerically these are about 0.02 and 0.15 
respectively a.t f3 — 6.0, and 0.01 and 0.1 respec- 
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tively at /3 = 6.2. 

The differences between our values and those 
determined by the ALPHA collaboration are con- 
sistent with this qualitative expectation. In fact, 
the differences between the two determinations of 
are 0.017(1) and 0.008(1) at /3 = 6.0 and 6.2 
respectively, exactly as expected from an O(a^) 
scaling. For Z% the differences of 0.020(17) and 
0.009(13) at the two /3 values are similar in mag- 
nitude, but have much larger errors. 

Of the coefficients bx, only by has been cal- 
culated by the ALPHA collaboration. The cor- 
responding differences, 0.14(5) and 0.13(2), have 
large errors but are not inconsistent with the ex- 
pected 0{a) scaling. 

The situation is less clear for the coefficients 
cx- The ALPHA collaboration has computed 
only ca and cy and our results are completely 
consistent with theirs at /3 = 6.2. The differences 
at (3 = 6.0 are much larger in comparison and 
marginally significant. 

5. Comparison with perturbation theory 

The perturbative results quoted in the table are 
based on one- loop tadpole improved perturbation 
theory and their errors are obtained by squaring 
the one-loop term in this tadpole improved se- 
ries. We notice that the non-perturbatively de- 
termined Zy and Zp/Zg are significantly lower 
than their perturbative values and all the 6x's 
are higher. The coefficients cx are small, and ex- 
cept for CA^ agree with perturbation theory within 
errors. 

It is worth mentioning that the difference be- 
tween the perturbative and non-perturbative Zy 
is about 0.046(1) at /? = 6.0 and 0.037(0) at 
/3 = 6.2, where the errors are purely statisti- 
cal. As an O(a^) contamination in our non- 
perturbative determination is expected to change 
by a factor of two between these two calculations, 
it cannot be a dominant contribution to the ob- 
served difference. Furthermore, noting that the 
tadpole improved ols is about 0.13 and 0.12 in 
the two calculations, these differences are only 
about 2.5 times a^, not unreasonable in a slowly 
converging perturbative series. The pattern is 
similar for Zp/Zg, where the difference between 



the perturbative and non-perturbative results is 
about 0.09(2) at /? = 6.0 and 0.07(2) at ^ = 6.2. 
Because of the larger errors, it is, however, not 
possible to rule out O(a^) artifacts in this case. 

Except for fey, the non-perturbative values for 
bx are constant within errors as we move from 
/3 — 6.0 to /3 = 6.2. The leading non-perturbative 
errors in these quantities are 0(a), which should 
change by a factor of about 1.4 between these two 
couplings. On the other hand, the perturbative 
one-loop term accounts for only about a fourth of 
by ~ \ and &^ — 1. Such a high estimate of bx with 
a correspondingly low value of Z\ could indicate 
a problem with understanding the quark mass de- 
pendences in general: more study is needed to 
clarify this issue. 

6. Conclusion 

We have demonstrated the feasibility of our 
method for determining the scheme-independent 
renormalization constants of the quark bilinear 
operators. Differences from previous calculations 
are consistent with the expected 0(a}) differ- 
ences. Perturbation theory, even when tadpole 
improved, seems to have 0(few x a^) residual er- 
rors for the chirally extrapolated renormalization 
constants. The improvement constants cx are 
small, but possibly somewhat larger than pre- 
dicted by perturbation theory; a detailed com- 
parison is still not possible due to the large er- 
rors. The mass dependences of the renormaliza- 
tion constants, given by bx^ are much larger than 
perturbation theory predicts. 
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